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1 INTRODUCTION 

Clusters of galaxies occupy an eminent position in the 
structure hierarchy, being the most massive virialized 
systems known and therefore they appear to be ideal 
tools for testing theories of structure formation and 
extracting cosmological informations (cf. Bahcall 1988; 
West Jones & Forman 1995; Bohringer 1995; Carlberg 
et al. 1996; Borgani & Guzzo 2001; Nichol 2002 and 
references therein). 

To investigate the global properties of the cosmo- 
logical background it is necessary to construct and study 
large samples of clusters (cf. Borgani & Guzzo 2001). 
This understanding has initiated a number of studies 
aiming to compile unbiased cluster samples to high red- 
shifts, utilizing multiwavelength data (e.g. optical, X- 
ray, radio). On the other hand, the study of individual 
clusters, provide complementary information regarding 
their physical properties and evolutionary processes. 
Overall, it is very important to fully understand the 
selection effects that enter in the construction of cluster 
samples since these could bias any statistical analysis of 
these samples (cf. Sutherland 1988). 

At optical wavelengths there are several available 
samples in the literature. For example, the Abell/ACO 
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catalogue (Abell, Corwin & Olowin 1989) was con- 
structed by visual inspection of the Palomar Observa- 
tory Sky Survey plates and is still playing an important 
role in astronomical research. Since then, a large num- 
ber of optically selected samples constructed with au- 
tomated methods have been constructed: EDCC (Edin- 
burgh Durham Cluster Catalogue; Lumsden et al. 1992) , 
APM (Automatic Plate Measuring; Dalton et al. 1994), 
PSCS (Palomar Distant Cluster Survey; Postman et 
al. 1996), EIS (ESO Imaging Cluster Survey; Olsen et 
al. 1999), RCS Red-Sequence Cluster Survey; Gladders 
& Yee 2000) and the Sloan Digital Sky cluster survey 
(Goto et al. 2002; Bahcall et al. 2003). The above clus- 
ter samples, based on different selection methods, aim 
to obtain homogeneously selected optical cluster sam- 
ples with redshifts that extended beyond the z ~ 0.2 
limit of the Abell/ACO catalogue. We should mention 
that the advantage of using optical data is the shear size 
of the available cluster catalogues and thus the statisti- 
cal significance of the emanating results. 

A major problem here is that the optical surveys 
suffer from severe systematic biases which are due to 
projection effects. Background and foreground galaxies, 
projected on the cluster could distort the identifications 
(e.g. Frenk et al. 1990), which is particularly true for 
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poor systems at high redshifts. X-ray surveys provide an 
alternative method for compiling cluster samples, owing 
to the fact that the diffuse Intra-Cluster Medium (ICM) 
emits strongly at X-ray wavelengths. This emission is 
proportional to the square of the hot gas density, result- 
ing in a high contrast with respect to the unresolved X- 
ray background, and thus X-ray selected clusters are less 
susceptible to projection effects. Therefore the main ad- 
vantage, inherent in the X-ray selection of flux-limited 
samples is that the survey volume and hence number 
densities, luminosity and mass functions can be reliably 
computed. Furthermore, X-ray cluster surveys can be 
used to study cluster dynamics and morphologies, the 
Sunyaev Zeldovich effect and finally their cosmological 
evolution. 

The first such sample, with large impact to the 
studies of clusters, was based on the Extended Ein- 
stein Medium Sensitivity Survey, containing 99 clusters 
(Stocke et al. 1991). The ROSAT satellite with its large 
field of view (FOV) and better sensitivity, allowed a leap 
forward in the X-ray cluster astronomy, producing large 
samples of both nearby and distant clusters (e.g. Cas- 
tander et al. 1995; Ebeling et al. 1996a, 1996b; Scharf 
et al. 1997; Ebelling et al. 2000; Bohringer et al. 2001; 
Gioia et al. 2001; Bohringer et al. 2002; Rosati, Borgani 
& Norman 2002 and references therein). Recently, the 
XMM-Newton observatory with ~ 10 times more effec- 
tive area and ~ 5 times better spatial resolution than 
the ROSAT provides an ideal platform for the study of 
galaxy clusters. 

However, even with the improved sensitivity of the 
XMM-TVewfon, optical surveys still remain significantly 
more efficient and less expensive in telescope time for 
compiling cluster samples, albeit with the previously 
discussed limitations (eg. incompleteness, projection ef- 
fects etc). Therefore it is necessary to study the differ- 
ent selection effects and biases that enter in detecting 
clusters in the two wavelength regimes. Donahue et al. 
(2002) using the ROSAT Optical X-ray Survey (ROXS), 
found that using both X-ray and optical methods to 
identify clusters of galaxies, the overlap was poor. About 
20% of the optically detected clusters were found in X- 
rays while 60% of the X-ray clusters were identified also 
in the optical sample. Furthermore, not all of their X- 
ray detected clusters had a prominent red-sequence, a 
fact that could introduce a bias in constructing cluster 
samples using only colour information (as in Goto et al. 
2002). 

The aim of this work is along the same lines, at- 
tempting to make a detailed comparison of optical and 
X-ray cluster identification methods in order to quan- 
tify the selection biases introduced by these different 
techniques. We use XMM-Newton which has a factor of 
~ 5 better spatial resolution and an order of magnitude 
more effective area at 1 keV, making it an ideal instru- 
ment for the detection of relatively distant clusters. 

The plan of the paper is as follows. The observed 
data sets are presented in Section 2. In Section 3 we 
discuss the methods employed to identify candidate op- 
tical clusters and comment on the systematic effects in- 
troduced in our analysis. Also, Section 3 describes our 
projected cluster shape determination method as well 



as the cluster surface brightness based on a King like 
profiles. In Section 4 we compare the optical and X- 
ray selected cluster samples. Finally, in Section 5, we 
present our conclusions. Throughout this paper we use 
Ho = 100ft km s^Mpc -1 and Q m = 1 - fi A = 0.3. 



2 OBSERVATIONS 

2.1 The optical data 

In our analysis we use the SDSS Early Data Re- 
lease (EDR), covering an area of ~ 400deg 2 in the 
sky (Stoughton et al. 2002). The SDSS is an ongoing 
imaging and spectroscopic survey that covers ~10000 
deg 2 of the sky. Photometry is obtained in five bands 
(u, g, r, i, z) to the limiting magnituge r < 22.5 mag, 
providing an homogeneous multi-color photometric cat- 
alogue. 

Goto et al. (2002), applied an object cluster find- 
ing algorithm to the SDSS EDR and produced a list 
of galaxy clusters (hereafter CE catalogue), with esti- 
mated photometric redshifts which contains 2770 and 
1868 galaxy clusters in the North and South slices, re- 
spectively. The cluster redshifts were estimated using 
color information by identifying the bin in g — r which 
has the largest number of galaxies around the color 
prediction of elliptical galaxies (Fukugita, Shimasaku, 
Ichikawa 1995) at different redshifts (which define the 
different g — r bins) . Note that the true and estimated 
redshifts are better correlated for z < 0.3, with rms 
scatter ~ 0.015, while for z > 0.3 it is ~ 0.021. Note, 
that the CE method is optimized up to z < 0.4 and 
becomes highly insensitive at higher redshifts. 

2.2 The XMM-Newton data 

We analyzed 9 XMM-Newton fields with nominal expo- 
sure times between 2 and 10 ksec. The XMM-Newton 
field-of-view is a circle with a radius of 15 arcmin and 
thus our survey covers an area of 1.8 deg 2 . How- 
ever, one of the fields, suffering from significantly el- 
evated and flaring particle background, was excluded 
from the X-ray analysis. This reduces our effective area 
to 1.6 deg 2 . The details of the X-ray observations are 
given in Georgakakis et al. (2004) and Georgantopoulos 
et al. in preparation. 



3 OPTICAL CLUSTER DETECTION 

Identifying real clusters using imaging data is a difficult 
task since projection effects can significantly affect the 
visual appearance of clusters. Many different algorithms 
have been proposed and applied on different data sets. 
In what follows we present a new cluster finding algo- 
rithm (based on two distinct algorithms) and apply it on 
the subsample of the r-band SDSS data*, that overlaps 
with our XMM- Netwon survey. 

* We have tested that our results remain qualitatively unal- 
tered when using the i-band. 
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Figure 1. The SDSS probability density function (pdf). 



3.1 The Smoothing+Percolation Procedure 

This cluster detection algorithm (hereafter SMP 
method) is based on smoothing the discrete distribu- 
tion using a Gaussian smoothing kernel on a N gr x N gr 
grid: 



W(| Xi -x gr |) = 



exp 



■2RL 



(1) 



with _R sm the smoothing radius in grid-cell units and Xj 
the Cartesian position of the i th galaxy. The smoothed 
surface density, of the j th cell at the grid-cell positions 



Pi ( x gr) — 



X; i Pj(x i )W(|x i -Xg 



/W(|x,- 



|)d 2 a: 
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where the sum is over the distribution of galaxies with 
positions Xi. For each cell we define its density fluctua- 
tion, 5j, as: 



8i = 



Pi( x gr) - (p) 

(p) 



(3) 



where (p) is the mean projected SDSS galaxy density 
and the probability density function, f(S), which is plot- 
ted in Fig. 1. We select all grid-cells with overdensities 
above a chosen critical threshold (S > S CI ) and then 
we use a friends-of- friends algorithm to link all adjacent 
cells in order to form groups of connected cells, which 
we consider as our candidate clusters. Note that the grid 
cell size is such that at z = 0.4 it corresponds to 100 
h^kpc (~ 19"). 

There are three free parameters in the above proce- 
dure: the critical overdensity threshold (<5 cr ) the smooth- 
ing radius (Rum) an d the number of connected cells (n c ) 
above which we consider a candidate cluster. Due to the 
fact that clusters should be identified as high density re- 
gions, we choose the critical value of the overdensity to 
be Scr — 1 which corresponds to the ~ 98 percentile of 



Figure 2. Statistical significance of our cluster detec- 
tion procedure as a function of the smoothing radius R Em 
and minimum number of connected cells, n c . The darkest 
greyscale corresponds to V = 0.96, while the faintest to 
V <0.88 . The continuous line corresponds to the expected 
number of APM clusters, while the broken line to the num- 
ber Goto et al. (2002) clusters. The cross reflects our best 
choice of the free parameters (-R S m and n c ). 



the probability density function (pdf) plotted in Fig.l. 
Had we used lower values of S CI we would have perco- 
lated overdense regions through large parts of the whole 
area. 

The Ram and n c parameters were chosen in such a 
way as to minimize the number of spuriously detected 
clusters. To this end we perform a series of Monte-Carlo 
simulations in which we randomize the positions of the 
SDSS galaxies, destroying its intrinsic clustering, while 
keeping unchanged the galaxy redshift selection func- 
tion. On this intrinsically random galaxy distribution 
we apply the procedure described above, by varying the 
values of the free parameters (f? sm and n c ). We define 
the probability of detecting real clusters in the SDSS 
data, as V(Ram, n c ) = 1— iV raa d/iVsDSS, were AT ran is the 
number of clusters detected in the randomized distribu- 
tion and iVsDSS is the number of clusters detected in the 
true SDSS data. In Fig. 2 we show as a greyscale these 
probabilities. Although, increasing i? sm and n c , results 
in very high probabilities of detecting real high-density 
regions, it also results in small numbers of detected clus- 
ters. We therefore attempt to break this degeneracy by 
using the expected number of clusters from existing sur- 
veys. For example, the continuous line in Fig. 2 corre- 
sponds to the expected number of APM clusters in the 
area covered by our XMM/2dF survey^, while the thick 



t This number (~ 40) is estimated, without assuming any 
evolution, by multiplying the local (z < 0.1) mean number 
density of APM clusters with the volume element covered 
by our survey out to a z ~ 0.45 within the concordance 
cosmological model (Q^ = 0.7). 
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Figure 3. The smooth SDSS density field on equatorial co- 
ordinates. The 5 C r = 1 level appears as a thick continuous 
line. The large dotted circles represents the XMM 15 arcmin 
radius fields of view of our shallow XMM/2dF survey. 



broken line delineates the corresponding number of the 
Goto et al. (2002) clusters. The region of the highest 
probability, that falls within the above cluster number 
limits, is that corresponding to: 



Rsm — 28.5 



n c = 12 , 



which we choose as our optimal parameters for the de- 
tection of real clusters. 

In Fig. 3 we plot the smoothed SDSS galaxy distri- 
bution, using the above value of i? sm , on the equatorial 
plane with a contour step of AS — 0.2, starting from 
S CI — 1. Many overdense regions are apparent. Finally, 
we note that our cluster finding algorithm with n c = 12 
produces a total of 25 candidate clusters, 11 of which 
are common with the Goto et al. (2002) sample of 29 
clusters in the same region. 



3.2 The Matched Filter Algorithm 

In addition to the SMP technique we have also em- 
ployed the matched filter algorithm (hereafter MFA) 
described by Postman et al. (1996) to identify optical 
galaxy overdensities. The advantage of this method is 
that it exploits both positional and photometric infor- 
mation producing galaxy density maps where spurious 
galaxy fluctuations are suppressed. Also, an attractive 
feature of this method is that it provides redshift es- 
timates for the detected candidate galaxy clusters. A 
drawback of the matched filter method is that one must 
assume a form for the cluster luminosity function and 
radial profile. Clusters with the same richness but differ- 
ent intrinsic shape or different luminosity function from 
the assumed ones do not have the same likelihood of 
being detected. 

A detailed description of the matched filter algo- 
rithm can be found in Postman et al. (1996). In brief, the 



filter used to convolve the galaxy catalogue is derived 
from an approximate maximum likelihood estimator ob- 
tained from a model of the spatial and luminosity distri- 
bution of galaxies within a cluster. At any given patch 
of the sky the number density of galaxies per magnitude 
bin is 



D(r, m) = b(m) + A c ; P(r/r c ) <j>(m — m*) 



(4) 



where D(r, m) is the galaxy surface density (clus- 
ter+field galaxies) at a given magnitude m and at a 
distance r from the galaxy cluster center, b(m) is the 
field galaxy surface density, P(r/r c ) is the cluster pro- 
jected radial profile, 4>{m — m*) is the cluster luminosity 
function and A c ; is an estimator of the cluster richness. 
The parameters r c and m* are the characteristic cluster 
scale length and apparent magnitude corresponding to 
the characteristic luminosity of the cluster luminosity 
function. Postman et al. (1996) showed that the filter 
for cluster detection should maximize the expression 

4>(m — 



I 



P(r/r c ) ■ 



b{m) 



D(m, r) d r dm 



(5) 



Assuming that the galaxy distribution D(m, r) can be 
represented by a series of (5-functions at the observed 
positions and magnitudes, the integral above is reduced 
to 



S(U) = £ P(r k /r c )L(m k ) , 



where L(mk) is defined as 



L(m) = 



<f>(m — m*) 
6(mj 



(6) 



(7) 



The quantity S(i,j) is the value of the (i, j) pixel of the 
convolved galaxy density map and Nt is total number 
of galaxies in the catalogue. L(mk) is the luminosity 
weighting function (i.e. flux filter). In practice, the sur- 
vey area is binned into pixels of given size (for 
the choice of values see below) and the sum in equa- 
tion 6 is evaluated by iterating through all galaxies in 
the catalogue. This then is repeated for every pixel of 
the density map. Both m* and r c are a function of red- 
shift and hence S(i,j) also depends on redshift through 
these parameters. The redshift dependence of m* also 
includes a fc-correction as discussed in the next section. 

The flux filter in equation 7 has a divergent integral 
at faint magnitudes for Schechter luminosity functions 
with a < —1. To overcome this problem Postman et al. 
(1996) modified the flux filter by introducing a power- 
law cutoff of the form 10" ' 4 Therefore, the flux 
filter in equation 6 is defined as: 



L(m) 



4>{m-m*)lQ- {)A ^ 
b(m) 



(8) 



The radial filter in equation 6 is defined: 



P(r)= — — 1 ; 1 , ifr<r co 

otherwise, (9) 

where r c is the cluster core radius and r co is an arbitrary 
cutoff radius. Here, we assume r c — 100fe _1 kpc and 
r co = 1 h^ 1 Mpc. The characteristic absolute magnitude 
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of the luminosity function, in r— band is taken to be 
M* = -20.83 + 51og/i (Blanton et al. 2003). Both the 
radial and the flux filter are normalized as described in 
Postman et al. (1996). 



3.2.1 The MFA Cluster detection 

The matched filter algorithm is applied to a 1.6 deg 2 
subregion of the SDSS centered on our XMM/2dF sur- 
vey. We consider galaxies with magnitudes brighter 
than r = 22.5 mag. At fainter magnitudes the SDSS 
is affected by incompleteness. The galaxy density map 
S(i,j) representing the galaxy density map is indepen- 
dently estimated for redshifts between z min — 0.1 to 
z m ax = 0.6, incremented in steps of 0.1. The z max cor- 
responds to the redshift where m* becomes comparable 
to the limiting magnitude of the survey. The character- 
istic luminosity, L* , the faint end slope of the luminosity 
function, a, and the cluster core radius, r c , are assumed 
to remain constant with redshift. 

The conversion from luminosity to apparent mag- 
nitude requires an assumption to be made on the k- 
correction of the galaxies. Here we assume a non- 
evolving elliptical galaxy model obtained from the 
Bruzual & Chariot (1993) stellar population synthe- 
sis code (Pozzeti, Bruzual & Zamorani 1996). ft is im- 
portant to(hereafter SMP method) emphasize that the 
choice of fc-correction does not significantly affect the 
cluster detections, but has an impact on the redshift 
determination. 

The pixel size of the galaxy density maps at any 
redshift is taken to be ~ 19" corresponding to a pro- 
jected cluster core radius of r c = 100/i _1 kpc at the 
redshift z — 0.4 (see section 3.1). Galaxy density maps 
are created for redshifts 0.1-0.6 and stored in FITS im- 
ages. The peaks of the galaxy distribution are detected 
using SExtractor (Bertin & Arnouts 1996). The mean 
and the variance of the background are determined us- 
ing a global value in each likelihood map. The main in- 
put parameters are the detection threshold, Odct, given 
as a multiple of background variance, and the minimum 
number of pixels, -/V min for a peak to be extracted as can- 
didate cluster. Simulations have been carried out (see 
next section) to optimize these parameters and to mini- 
mize the number of spurious cluster detections (« 5%). 
We adopt a dct = 4.0 while iV m i n is set to the area of a 
circle with radius r c — 100 h~ x kpc at any given redshift. 

The significance of a detection, a, is defined as the 
maximum value of all the pixels associated with that 
detection (i.e. pixels with values above the SExtractor 
detection threshold), expressed in multiples of the rms 
noise above the background. The peaks (i.e. candidate 
clusters) extracted on density maps of different redshifts 
are matched by positional coincidence using a radius of 
~ 1 arcmin. The approximate redshift of the candidate 
cluster is taken to be the redshift where the detection 
significance is maximum. 

Within the area of our XMM/2dF survey the MFA 
method identifies 12 cluster candidates out of which 6 
are found also by Goto et al. (2002) and 9 are found by 
the SMP method. 



3.2.2 The MFA simulations 

Monte Carlo simulations similar to those described in 
section 3.2 were employed to test the performance of 
the matched filter algorithm and to establish the best 
choice of extraction parameters. A total of 100 random 
galaxy catalogues over an area of ~ 1.6 deg 2 (i.e. similar 
to that covered by the real catalogue) were produced 
with the same magnitude distribution as the real data- 
set. The matched filter algorithm was then applied to 
these mock catalogues and the SExtractor was used to 
identify peaks in the derived density maps. 

Following Olsen et al. (1999) we minimize the num- 
ber of spurious detections using the SExtractor param- 
eters iV m i n and adct as well as other properties of the 
noise peaks such as the number of redshifts where a 
peak is detected, n z . We adopt <7d c t = 4, iV m i n = 7rr 2 at 
any given redshift and we only consider peaks that ap- 
pear in at least two consecutive redshifts (i.e. n z > 2). 
For this choice of parameters we estimate the spurious 
cluster contamination to be w 5%. 



3.3 Our final optically selected cluster sample 

We construct our final optical cluster catalogue by 
adopting the conservative approach of considering as 
cluster candidates those identified by both the SMP 
and MFA techniques. Our final sample, which we call 
SMPMFA, consists of 9 clusters with SDSS richness of 
more than 20 galaxies, corresponding roughly to APM 
type clusters. This estimated richness corresponds to the 
number of galaxies having r-band magnitudes < m-j + 2 
(where m-i is the magnitude of the third brightest clus- 
ter member) within a radius of 0.75 h^ 1 Mpc, estimated 
using the approximate redshift of the cluster. 

Comparing our SMPMFA sample with the 29 clus- 
ters of Goto et al. (2002), we find 5 in common. 



3.3.1 Cluster Shapes 

We investigate the reality of our candidate clusters by 
fitting a King's profile to their projected galaxy distri- 
bution. We expect that a galaxy enhancement produced 
by projection effects should have a rather flat profile and 
thus the goodness of fit should be low^ . The King's pro- 
file is given by: 



E(0) tx 



1+ ^ 



(10) 



where 9 C is the angular cluster core radius. The slope 
a is related to /3, the ratio of the specific energy in 
galaxies to the specific thermal energy in the gas, by 
a — (3/3 — l)/2 and spans the range 0.6<J a<J 1 (cf. 
Bahcall & Lubin 1993; Girardi et al. 1995; Girardi et al. 
1998). 

In order to quantify the parameters of the fit we 
perform a standard x 2 minimization procedure between 



i Note however that a real cluster with significant substruc- 
tures could also produce a rather poor King's profile 
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Table 1. List of the SDSS clusters in the XMM area. The correspondence of the columns is as follows: index number, right 
ascension a and declination <5 of the cluster center, cluster ellipticity e and <f> is the cluster position angle (in degrees), the cross 
correlation results (C-C)- the numbers correspond to: 1 - CE method, Goto et al. (2002); 2 - Matched Filter Algorithm; 3 - 
X-ray extended sources (Gaga et al. in preparation); 4 - Couch et al. (1991). The next columns are: the estimated redshift by 
different methods, the core radius with its 2a error and the x 2 probabilities (P x 2) of consistency between the intrinsic cluster 
density profile and King model. Finally, Af g , 7713 and f x are the number of galaxies, the magnitude of the third- brightest cluster 
member (within a radius ~ 0.75 h~ 1 Mpc) and the detected flux or the 3a upper flux limits in units of 10 — 14 erg cm -2 s — 1 
respectively (for details see section 4). Note that the r c has units of h~ 1 Mpc and the King profile slope used a ~ 0.7. 



Index 


Q 


8 


e 




C-C 


z 


r c 


P X 2 


M g 


7713 


fx 


1 


205.513 


-0.316 


0.10 


113.7 


2 


0.10 2 


0.10 ±0.03 


0.14 


29 


17.38 


<0.40 


2 


205.417 


0.287 


0.56 


6.3 


1,2,3 


0.39 1 - 0.40 2 


0.50 ±0.17 


0.003 


26 


18.96 


1.86 


3 


205.307 


-0.183 


0.61 


4.7 


1,2 


0.40 1 - 0.60 2 


0.39 ±0.13 


0.73 


34 


19.37 


<0.16 


1 


205.550 


0.378 


0.37 


170.4 


1)2 


0.31 1 - 0.50 2 


0.52 ±0.17 


0.95 


19 


19.76 


<3.91 


5 


205.859 


0.354 


0.37 


90.7 


2 


0.60 2 


0.56 ±0.19 


0.50 


50 


20.06 


<6.37 


6 


205.777 


-0.022 


0.10 


158.2 


2,3,4 


0.60 2 - 0.67 4 


0.35 ±0.11 


0.98 


33 


20.03 


5.55 


7 


206.296 


0.078 


0.56 


22.0 


1,2 


0.37 1 - 0.30 2 


0.42 ±0.14 


0.59 


30 


19.25 


<7.29 


8 


206.363 


-0.438 


0.43 


3.3 


2 


0.30 2 


0.46 ±0.15 


0.63 


29 


19.49 


<3.29 


9 


206.210 


0.468 


0.10 


48.0 


1,2,3 


0.36 1 - 0.30 2 


0.09 ±0.03 


0.46 


28 


19.56 


6.03 



the measured surface density of each of our clusters and 
eq.(10): 



= E 



^data 



(11) 



In Table 1 we present our cluster positions and the pa- 
rameters of the King's profile, as well as the \ 2 probabil- 
ity (P x 2 ) of consistency between the King's profile and 
the cluster density profile. The errors are 2a estimates, 
corresponding to A\ 2 = 6.17 (A\ 2 = X 2 — Xmin, with 
Xmin the absolute minimum value of the \ 2 )- All but 
one candidate clusters are well fitted by a King's profile 
with P x 2 > 0.14, while the slope a of the profile is ~ 0.7, 
corresponding to a /3-parameter of ~ 0.8 (cf. Cavaliere 
& Fusco-Femiano 1976; Jones & Forman 1999). Further- 
more, we present the estimated cluster redshifts (either 
from Goto et al. 2000 or from the MFA method). In the 
case of cluster # 6 we give the spectroscopic redshift 
from Couch et al. (1991). 

For those clusters which are well fitted by a King- 
like profile we can give an approximation regarding their 
spatial core radius (r c ). Different studies (Girardi et al. 
1995; Girardi et al. 1998), have shown that the radial 
core radius of Abell clusters has an average value of 
f c ~ 0.10 — 0.15 ft -1 Mpc. Therefore, having an estimate 
of the cluster redshift, the corresponding core radius can 
be easily found from r c = dLtan^c), where g!l is the 
luminosity cluster distance. 

In order to check for possible evolutionary trends 
regarding the cluster core, we have plotted the core ra- 
dius as a function of redshift. As shown in Fig. 4 we find 
a significant correlation; the coefficient is ~ 0.60 with 
the probability of zero correlation being < 10 . In the 
local universe (0 < z < 0.15), the core radius is approx- 
imately equal with ~ 0.10ft _1 Mpc, which is similar to 
the value derived by Girardi et al. (1998). The line of 
the best fit is given by 



0.530(±0.291)2 + 0.138(±0.134) 



(12) 




In order to check if the observed evolutionary trend 
is an artifact of possible systematic effects, we have 



Figure 4. The cluster core radius (r c ) redshift correlation 
for the analyzed SMPMFA. In the insert we show the SDSS 
redshift selection function using the Blanton et al. (2003) 
luminosity function. 



firstly investigated the effect of the SDSS redshift se- 
lection function. We have determined the predicted 
SDSS redshift distribution using the r-band luminosity 
function of Blanton et al. (2003) for galaxies between 
0.02 < z < 0.22. We have extrapolated their luminos- 
ity dependent density evolution model up to z = 2 and 
find that the SDSS redshift distribution increases up to 
z ~ 0.45 (insert of Fig. 4). This implies that the sam- 
ple is roughly volume limited out to this distance and 
thus the r c — z trend seen in Fig. 4 should not be due to 
systematic effects related to undersampling of the clus- 
ter galaxy population. Since we have used a single band 
(r-band) to detect clusters, we also explore the effect of 
determining their morphological characteristics in dif- 
ferent rest-frame bands. We note that the wavelength 
difference (effective A) between the Sloan r and z bands 
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is ~2800 A, which corresponds to a redshift difference of 
~0.45. Therefore, the rest-frame r-band at low redshifts 
(z < 0.15) roughly shifts to the z-band at z ~ 0.6. We 
then compare the mean core-radius of our most distant 
clusters (at ~0.6 and 0.67) using the z-band (within 
its completness limit m z < 21.5) with the mean core- 
radius of our relatively nearby clusters (z < 0.3) deter- 
mined in the r-band (within its the completness limit of 
m r < 22). We find that the r-band core radius of the 4 
low-z clusters is 0.27 ±0.19 while the z-band core radius 
of the 2 high-z clusters is 0.65 ±0.35. Therefore it seems 
that the effect of the different rest-frame bands in the 
definition of the cluster morphology (at least for our 
sample and the redshifts covered) does neither create 
nor mask the observed evolutionary trend. 

We list in Table 1 the cluster shape parameters 
determined by using the moments of inertia method 
(cf. Carter & Metcalfe 1980; Basilakos, Plionis, Mad- 
dox, 2001). This method is based in evaluating the mo- 
ments In = E <5j(r] - a%), I22 = E Sj(r] - y]), 
I12 = I21 — — E &i x iVu where 8j is the galaxy overden- 
sity of each cell and x, y are the Cartesian coordinates of 
each cell, after transforming the equatorial coordinates 
into an equal area coordinate system, centered on the 
cluster center. Then, diagonalizing the inertia tensor 

det(Iy - A 2 M 2 ) = (M 2 is 2 x 2 unit matrix.) (13) 

we obtain the eigenvalues Ai, A2, from which we define 
the ellipticity of the configurations under study by: e = 
1 — A2/A1, with Ai > A2. The corresponding eigenvectors 
provide the direction of the principal axis. It is evident 
that 3 out of 9 clusters have large projected ellipticities 
(e > 0.5). 



4 COMPARING CLUSTERS IN THE 
OPTICAL & X-RAY BANDS 

4.1 X-ray Cluster Detection 

The full details of the cluster detection procedure and 
the cluster X-ray luminosity and temperature determi- 
nations will be described in Gaga et al. (in preparation). 
Here we only briefly present the main method and re- 
sults. 

In order to detect candidate clusters in our XMM 
fields we use the soft 0.3-2 keV data since this band 
maximizes the signal to noise ratio especially in the 
case of relatively low temperature galaxy clusters. In 
particular we use the EWAVELET detection algorithm 
of the XMM-Newton SAS v. 5. 3 analysis software pack- 
age, which detects sources on the wavelet transformed 
images. 

We search for sources down to the 5a detection 
threshold, on the PN and the co-added MOS1 and 
MOS2 detector images, separately. The output list is fed 
into the SAS EMLDETECT algorithm which performs 
a maximum likelihood PSF fit on each source yielding a 
likelihood for the extension. We are using an extension 
likelihood threshold corresponding to a level of less than 
~ 0.5% spurious extended sources. We have detected 7 
candidate clusters on the MOS mosaic, while 5 extended 
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Figure 5. The candidate cluster positions within our shallow 
XMM- Newton survey. Open squares are the optically selected 
clusters using the SMPMFA technique while open circles are 
the Goto et al. (2002) clusters. Crosses are X-ray selected 
clusters of Gaga et al. (in preparation 2004). The large dotted 
circles represent the 15 arcmin radius XMM fields of view 
covered by our shallow XMM survey. 

sources were detected on the PN images, out of which 
3 overlap with the MOS candidates. 

We have also visually inspected the candidate clus- 
ters and found that in some occasions the high extension 
likelihood was due to multiple point sources (confusion). 
These spurious detections are excluded from our final X- 
ray selected cluster candidate list. This results in a final 
list of 4 X-ray cluster candidates, out of which 2 are de- 
tected on both the MOS and PN, one is detected only 
on the MOS (there was no valid PN image) and the last 
one is detected only on the PN (note that this cluster 
was detected with a slightly smaller likelihood probabil- 
ity, ~ 0.99). The faintest extended source has a flux of 
~ 2 x 10" 14 erg cm" 2 s" 1 . Note that for the (0.3-2) keV 
energy band the net X-ray count rate for our detected 
extented sources is 3 to 7 x 10 _3 cts s _1 , within a 30 
radius circle. 

In Table 2 we present the properties of the X-ray 
cluster candidates. 

4.2 The optical versus X-ray cluster 
detections 

We perform a cross-correlation between the different 
cluster catalogues (optical, X-rays), using a 3 arcmin 
matching radius to identify common objects. In Fig. 
5, we plot the positions of (a) our optically detected 
cluster candidates (SMPMFA sample), (b) the Goto et 
al. (2002) clusters and (c) the X-ray detected XMM- 
Newton clusters. All four X-ray detections coincide with 
optical cluster candidates from different methods, with 
the largest coincidence rate (3 out of 9) being with our 
SMPMFA optical candidates. The most distant cluster 
in our SMPMFA sample, located at a =13h 43min 5 sec 
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Table 2. The X-ray cluster properties: (1) and (2) equitorial coordinates (J2000) (2) X-ray Luminosity (erg s -1 ) in the 0.3-2 
keV band, (3) X-ray flux (0.3-2) keV, in units of 10 _14 erg sec -1 , cm -2 , (4) cluster temperature from Gaga et al. (in preparation 
2004), (5) the XMM field, (6) optical follow up and finally, the cluster rcdshift. Note that the subscripts 1,2 and 4 are the same 
as in Table 1. 



a 


S 




h 


T(keV) 


XMM-Field 


Optical 


z 


detector 


13 41 39.1 


+00 17 39.3 


10 42.63 


1.86 


~3 


F864-1 


SMPMFA+CE 


0.40 2 


PN+MOS 


13 43 04.8 


-00 00 56.3 


I p43.60 


5.55 


~ 2 


F864-5 


SMPMFA 


0.67 4 


PN+MOS 


13 45 11.9 


-00 09 52.6 


10 42.00 


5.05 


~ 1 


F864-6 


CE 


0.12 1 


MOS 


13 44 49.7 


+00 27 53.3 


10 42.87 


6.03 


~ 1 


F864-3 


SMPMFA+CE 


0.30 2 


PN 



Table 3. Cross correlations results. The correspondence of the columns is as follows: comparison pair and the number of 
common clusters. 



Methods compared 


Number of common clusters 


SMPMFA - CE 


5 


SMPMFA - X-rays 


3 


CE - X-rays 


3 



S = 00° 00 56 , found also in X-rays, is missed by Goto 
et al. (2002). 

The cross-correlation results between the different 
cluster detection methods are presented in Table 3, 
which in the first column lists the comparison method 
pair and in the second column the number of common 
objects. 

It is interesting that our SMPMFA combined 
method finds 3 out of the 4 X-ray clusters, missing the 
nearest one, which from its X-ray luminosity appears 
to be a group of galaxies, rather than a cluster. How- 
ever, there are still 6 SMPMFA clusters that do not 
appear to have X-ray counterparts, for which we give 
their 3a flux upper limits (see Table 1 last column). 
Conversion from count rate to flux was done assuming 
a bremsstrahlumg spectrum with 3 keV. The flux was 
measured in a 30 radius cell and the correction to the 
total flux was performed taking into acount a King pro- 
file with r c = 150 h~ 1 kpc. This could be a hint that 
these clusters are either the results of projection effects 
(which however cannot explain the good King profile 
fits - see Table 1), or that our XMM-Newton survey is 
too shallow to reveal the probably weak X-ray emission 
from these clusters. 

In order to address this issue and to study the rela- 
tion between the limiting flux of our X-ray survey with 
respect to exposure time, we have carried out the follow- 
ing experiment. We have analysed observations taken 
from 15 XMM-Newton public fields with mean expo- 
sure times ~ 21 ksec in the soft 0.3-2 keV band, after 
filtering to correct for high particle background periods. 
Using the parameters of the SAS software as described 
previously, we have detected 31 candidate clusters. The 
faintest cluster detected, with a flux of ~ 5 x 10 -15 erg 
cm~ 2 s^ 1 , was found in the deepest field with an expo- 
sure time of 37 ksec. We then reduce the exposure times 
to a new mean value of ~ 5 ksec, similar to our shallow 
survey, and find only 9 out of the 31 previously iden- 
tified candidate clusters (29%), having a limiting flux 
of ~ 2.3 x 10~ 14 erg cm~ 2 s _1 . Therefore, had we had 



deeper XMM-Newton observations (by an average factor 
of ~ 5 in exposure time) we would have detected ~ 13 
X-ray candidate clusters in the region covered by our 
shallow XMM-Newton/2dF survey, which is consistent 
(within la) with the number of our SMPMFA cluster 
candidates. 



5 CONCLUSIONS 

We have made a direct comparison between optical 
and X-ray based techniques used to identify clusters. 
We have searched for extended emission in our shallow 
XMM-Newton Survey, which covers a ~ 1.6 deg 2 area 
(8 out of 9 original XMM pointings) in the North Galac- 
tic Pole region and we have detected 4 candidate X-ray 
clusters. We have then applied a new cluster finding 
algorithm on the SDSS galaxy distribution in this re- 
gion, which is based on merging two independent selec- 
tion methods - a smoothing-percolation SMP technique, 
and a Matched Filter Algorithm (MFA). Our final opti- 
cal cluster catalogue, called the SMPMFA list, counts 9 
candidate clusters with richness of more than 20 galax- 
ies, corresponding roughly to the APM richness limit. 

Out of the 4 X-ray candidate clusters 3 are common 
with our SMPMFA list. This relatively, small number 
of optical SMPMFA cluster candidates observed in X- 
rays suggest that some of the optical cluster candidates 
are either projection effects or poor X-ray emitters and 
hence they are fainter in X-rays than the limit of our 
shallow survey f x (0.3-2keV) ~ 2 x 10~ 14 erg cm" 2 s" 1 . 
This latter explanation seems to be supported from an 
analysis of public XMM-Newton fields with larger expo- 
sure times. 
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